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Abstract 

An in-depth theoretical study is carried out to examine the quasi-deterministic entanglement of 
two atoms inside a leaky cavity. Two A-type three-level atoms, initially in their ground states, 
may become maximally entangled through the interaction with a single photon. By working out an 
exact analytic solution, we show that the probability of success depends crucially on the spectral 
function of the injected photon. With a cavity photon, one can generate a maximally entangled 
state with a certain probability that is always less than 50%. However, for an injected photon with 
a narrower spectral width, this probability can be significantly increased. In particular, we discover 
situations in which entanglement can be achieved in a single trial with an almost unit probability. 



PACS numbers: 03.67.Mn, 03.65.Ud, 42.50.Ct 



I. INTRODUCTION 



Entanglement has long been considered a fundamental property in quantum the- 



ory 



It has gained renewed interest due to its potential applications in quantum 



computations and information in recent years, and is generally considered as the corner- 



stone of this field (see, e.g., Refs. and references therein). While entanglement 

can be generated by a lot of different implements such as 
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: entanglement 
12 1, and ion trap jl3L Il4l 15 1, it i 



If 



is most 
However, 



nuclear magnetic resonance in liquid samples 

conveniently realized by photons, e.g., in parametric down-conversion 
though being ideal carriers of quantum information, photons cannot serve as convenient 
storage media. The storage of information is much better implemented with atoms. In this 
sense, cavity quantum electrodynamics (CQED) provides a perfect stage for the interaction 



of information carriers and depository to take place 
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2JJ. With the advance in CQED 



technologies, some proposals have also been made to generate entangled states of atoms in 



optical and microwave cavities 
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26] , many of which employ the leakage of 



the cavity not only as a channel for information passage, but also as the crucial element for 
^ entreat HQ. „ a ddMon , .ate paction tec^ using quantum 
measurement have also been proposed to generate entanglement in a probabilistic manner 



24. 



26]. 



Recently, Hong and Lee have proposed a scheme to generate entangled atoms inside a 



leaky cavity in a "quasi-deterministic" manner 



in Ref. 



2811 . In their model, which is similar to that 



2l| . two A-type three-level atoms, initially prepared in their ground states, may 



become maximally entangled after interacting with a single cavity photon. The polarization 
of the output photon, serving as an indicator of the atomic state, is then measured, and 
generation of entangled atomic state is successful if the photon carries the right polarization. 
Adopting the master equation approach, they numerically showed that under the optimal 
condition in which the ratio of the coupling strengths of the two channels is a/2, the maximum 
probability of success in each operation is 50%. They further argued that a unit probability 
of success can be achieved by introducing the output photon to the leaky cavity repeatedly. 
Through the application of this automatic feedback scheme, the two atoms in the cavity are 
entangled in a "quasi-deterministic" manner. 

The method proposed by Hong and Lee is appealing. However, their scheme is still 
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probabilistic for each single trial. It is therefore interesting to study whether the probability 
of success in each trial could be increased. The master equation approach adopted in Ref . 28| 
has limited further development in this direction. In addition, they obtained the optimal 
condition and hence the maximum probability by direct inspection on the numerical results. 
In the present paper, we first develop a rigorous analytic solution to their model and in turn 
provide a physical interpretation for the optimal condition mentioned above. We are then 
able to show explicitly that the upper bound 50% of the probability of success is an intrinsic 
nature of their scheme that makes use of a cavity photon. 

Furthermore, by generalizing the initial photon state to consider photons injected into 
the cavity with certain spectra, we find that the probability of success in each trial can be 
increased and larger than 50%. It is worthwhile to note that the master equation approach 
employed by Hong and Lee in Ref. 28( is not a proper description for the injected and 



output photons, which are usually characterized by wave packets with certain durations or, 
equivalently, by their respective spectral functions. The master equation approach, on the 
other hand, deals mainly with the "decay" of photons existing (or created) in the cavity. A 
photon that has already leaked out of the cavity is considered as a "decay" event and is no 
longer part of the system. As a consequence, their approach fails to describe the quantum 
states of the injected and output photons and is therefore deemed inappropriate in this 
situation. 

In this paper, we adopt the continuous frequency mode approach that considers the 
cavity and its environment as a single entity to quantize the entire system 2^], 0, l^ . 
This enables us to study a photon with an arbitrary spectral function. We find that the 
probability of success in each operation depends crucially on the spectral function of the 
injected photon. In general, the probability can be increased to well over 50%. As long as 
spontaneous decay can be neglected, with suitably chosen spectral function of the photon, 
deterministic entanglement generation can be realized irrespective of the leakage rate of 
the cavity. Remarkably, an almost unit probability of success can be obtained for realistic 
CQED setups and photon packets. 

The structure of our paper is as follows: In Sec. |Hj we introduce a proper set of spatial 
mode functions — the continuous frequency modes, in terms of which the electromagnetic 
field is quantized. The interaction Hamiltonian and the evolution of the atoms governed by 
it are discussed in Sec. IIHI and IIV| respectively. In Sec. El the generation of entanglement 
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is studied for two cases with different input states, namely, a cavity photon 32, |33j, |3J] and 
an arbitrary injected photon. We conclude the main features of our scheme in Sec. IVII 



II. CONTINUOUS FREQUENCY MODES 

In the following discussion we adopt the continuous frequency mode quantization scheme 
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321 ] . Hence 



to analyze optical processes taking place inside a leaky optical cavity [29, 
photons are represented by quantum states \ki) and |/cr), which are characterized by a 
continuous wave number k, and the polarization (L and R) of the photon. For defmiteness 
and without loss of generality, we consider in this paper a one-sided optical cavity as shown 
in Fig. which has length / and is bounded by two mirrors. The left mirror at x = 
is perfectly reflecting; the right mirror at x — I is partially transparent and characterized 
by frequency-dependent reflection and transmission coefficients r(k) and t(k). The spatial 



mode functions Uk(x) corresponding to quantum states |/cl(_r)) are given by 

I(k) sin kx < x < L 

U k (x) = { . 1 " (2.1J 

e -ikx _|_ ^{J^^x I < X < OO, 



where 



I(k) = r-^L, (2.2) 



1 + re 
-r — t + re 



-2ikl 



R ( k ) = x + re 2ua ■ ( 2 - 3 ) 

Notice that = 1, as required by conservation of energy. Standard Sturm-Liouville 

type orthogonality integral leads to the general normalization condition 

e(x)U k (x)U;,(x)dx = 2n5(k - k'), (2.4) 



o 

where e(x) is the position-dependent relative permittivity. For example, if an infinitely thin 
dielectric mirror is placed at x = I, then e(x) = 1 + (5(x — I), where ( determines the finesse 
of the mirror 

The quantization of the field is hence accomplished by decomposing the field operator 
A(x, t) into these continuous frequency modes and by introducing the usual annihilation 
(creation) operator (ojL), with fi = L,R being the polarization index. These operators 
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satisfy the standard commutation relation [d k ^, a\, ,] = 5^i5(k — k'). The result reads: 



A(x,t) 



fi=L,R 



k 1/2 \u k (x)a klx e 



- tkt + U* k (x)ale lkt 



dk. 



(2.5) 



Hereafter we adopt the cgs units and take h = c = 1. 

The coefficients I(k) and R(k) possess poles k n , which are roots of the equation \ J r re 2lkl = 
0. As long as r(k) and t(k) are slowly varying functions of k, the roots of the above equation 
are approximately given by: k n = k n — in n /2, where 

(2n + 1)tt -arg[r(fc n )] 



k r , 



Hi. 



21 

log |r(fc w )| 
I 



(2.6) 
(2.7) 



and n = 0, ±1, ±2, .... These complex frequencies define the quasi-modes of the cavity 
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31 



29 



35l | . Figure El shows a sketch of \U k {xo)\ 2 versus k at a point Xq inside the cavity. Each 
peak in the figure corresponds to a quasi-mode of the cavity. For large n, the separations 
and widths of the modes can be taken as constants given by tt/1 and k, respectively. The 
physical significance of these quasi-modes is that for a "good cavity", where |r| ~ 1, the 
mode function is negligibly small except at frequencies k ~ k n . In fact, it can be shown 



that the e! 
modes 



ectrodynamics inside the cavity is completely describable in terms of these quasi- 
. In the following, we shall make use of this formalism to discuss the interaction 



between the atoms and the quantized photon field. 



III. ATOM-FIELD INTERACTION 



We consider two identical A-type atoms A and B located near the center of a one-sided 
cavity as shown in Fig. ^ The separation of the atoms is assumed to be much smaller than 
the wavelength of the field in interest. The two ground states and the excited state are, 
respectively, denoted by \L a ), \R a ) and \e a ), where a — A, B refers to atom A and B. 
The excited state couples with the ground state \L a ) (\R a )) by emitting a \ki) (\k R )) mode 
photon. In our model, we assume that the ground states are stable, and the spontaneous 
decay rate of the excited states into non-cavity modes is negligible compared with the cavity 
leakage rate k. We shall use the notation \\xa\^b] k^) for states with two ground-state atoms 
and a single photon, where fiA, ^b, H = L, R and the first and second positions are assigned 



5 



to atom A and B, respectively. States with an excited-state and a ground-state atom will 
similarly be denoted by \eA^B\ 0) and \ha&b\ 0), where |0) is the vacuum state of the field. 

Assuming that the two ground states are degenerate and the energies of the ground and 
excited states are zero and u) e , respectively, in terms of the continuous frequency mode basis, 
the Hamiltonian can be written as 

H = UJ e \e a )(e a \+ dk Y ^a^a^ 

a=A,B J ° n=L,R 

poo 

+ dk Y 9 fl (k)ak^\e a ){ijL a \+h..c., (3.1) 

-'O „,_ A O 



a=A,B 
fj,=L,R 



where the function g{k) is proportional to the mode function Uk(x) evaluated at the position 
of the atoms, x a ~ 1/2. 

In the present paper, the atomic frequency is assumed to be close to one of the cavity 
quasi- mode frequencies given by Eq. (j2.6j) . Therefore, one can expand the mode functions 
around the resonance frequency and obtain the so-called single mode approximation result 



aiuium 

0Q, r 



where k c and k are, respectively, the frequency and the decay rate of the quasi-mode in 
resonance, and 6^ is a trivial phase angle. Besides, 



K= \9,{Wdk (3.3) 

J — oo 

is a measure of the coupling strength that depends on the dipole moment and the location 
of the atom. Also notice that within the single mode approximation, the lower limit of the 
frequency domain will hereafter be extended from to — oo, as shown in Eq. (|3.3|) . 



IV. ANALYTIC SOLUTION 

The system is initially prepared in the state 

/oo 
dk'f(k')\LL;k' L ), (4.1) 
-oo 

where f(k') is the spectral function of a one-photon state in the continuous frequency mode 
representation, satisfying the normalization condition: 

/oo 
\f(k')\ 2 dk' = l. (4.2) 
-oo 
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The explicit form of the injected photon will be specified later. 

We introduce an essential state basis to simplify the dynamics of the two atoms as follows. 
The frequency and polarization of the photon are two independent degrees of freedom of the 
Hilbert space. Restricted to single excitation, one of the two atoms may be excited while 
the other can be in \L) or \R), and thus there are four distinct states with a single excited 
atom. When the two atoms are in the ground states (\L) and \R) states), the polarization 
of the emitted photon can be either L or R, giving rise to a total of eight distinct states for 
each frequency. However, if we are only interested in the evolution with the specific initial 
state given by Eq. (|4.1j) . the dimension of the subspace involved is greatly reduced. One of 
the atoms in the initial state \LL; k' L ) may absorb the photon and the state evolves into 

\E) = ^=(\eL;0) + \Le;0)). (4.3) 

When the excited atom de-excites by emitting another photon of wave number k (k ^ k' in 
general), the resulting state may be \LL; k L ) or |$; k R ), where 

|$) = {\LR) + \RL)) (4.4) 

is a Bell state of the two atoms — the goal of the current discussion. Hence, the dynamics 
of our system is adequately described by transitions between the excited state \E) and the 
single-photon states \LL; ki) or |<&; kit). Considering hereafter an initial state specified by 
Eq. ([4.1)1 . we obtain an effective Hamiltonian within the above-mentioned subspace: 

/■oo 

H ea = u e \E}(E\ + / dkk(\LL;k L )(LL;k L \ + \$;k R )($;k R \) 

J —oo 

+ [°° dk (V2g L (k)\E)(LL; k L \ + g R (k)\E}($; k R \) + h.c. (4.5) 

J — oo v ' 

It is noteworthy that the dynamics of the current system, governed by H e s, bears strong 
resemblance to that of a single A-type three-level system in an optical cavity, with coupling 
strengths of the two channels being in the ratio of y/2X L to \ R . The factor y/2 can be 
explained as follows. First, while both atoms in state \LL; ki) can interact with the photon, 
one of the two atoms in state |$; k R ) is only a spectator. Second, the collective quantum 
effect of the state |$) contributes an enhancement factor y/2. These two effects together 
give an overall y/2 factor. Hence, if the coupling strengths and A# are in the ratio of 
1 : y/2, satisfying the optimal condition discovered in Ref . 2^ , the original two-atom model 
reduces to a symmetric three-level system. Interestingly enough, it is well known that a 



resonantly driven symmetric A-type atom located in an ideal cavity can exhibit perfect 
oscillations between its two ground states. We shall see that these analogies do shed light 
on our current investigation. 

To simplify the calculations, we perform another basis transformation to reduce the model 
to a two-level system. By defining 

\Uk)) = (V2g* L (k)\LL; k L ) + g R (k)\$; A*)) , (4.6) 

\Mk)) = ^(g R {k)\LL;k L )-V2g L {k)\$;k R )), (4.7) 

where 



V(k) = y/2\g L (k)\* + \g R (k)\*, (4.8) 

one of the ground states is turned into a dark state \ip2{k)), which does not take part in the 
interaction. The Hamiltonian can thus be written as 

H eS = H + V + H daxk , (4.9) 

where 

H = u e \E){E\ + / dkk\Mk))(Mk)\, (4.10) 



V= / dkV(k)\E)(ip 1 (k)\ + h.c, (4.11) 

J —oo 

#dark = / dkk\1> 2 (k))(Mk)\, (4-12) 

J — oo 

with -f^dark being the free Hamiltonian of the dark states \ip2{k))- 

To study the output state |^) ut of the system subject to an arbitrary state at t = 0, 
we adopt the resolvent method (see, e.g Ref. [3^) to deal with transitions between states 
with different /c's. Remarkably, if the trivial evolution of the dark state is discarded, the 
dynamics governed by the effective Hamiltonian is ostensibly analogous to that of a two-level 
atom interacting with a single quasi-mode. To simplify the calculation, the dark states will 
be ignored when solving the system using resolvent method. 

The resolvent of the Hamiltonian Hq + V is given by 

6{u>) = (4.13) 

with uj being a complex variable. It yields the retarded Green's function, 9(t) exp[— i{Hq + 
V)t], through an integral transformation: 

9(t)exp[-i(H + V)t] = lim — / G{u)e~ iwt du. (4.14) 
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In the following discussion, unless otherwise stated, we assume t > and therefore identify 
the retarded Green's function with the evolution operator U(t, 0) = exp[— i(Ho + V)t\. 

To evaluate the matrix elements of the resolvent, we partition the Hilbert space into two 
complementary parts with projection operators P = \E)(E\ and Q = 1 — P. The projection 
of G{uS) by P is given by 36]: 

P 



PG(uj)P 



UJ 



PH P - PR(uj)P 
where R(uj) is the level-shift operator defined as 

Q 



R( u ) = v + V 
It is then straightforward to show that 

(E\R(w)\E) 



uj - QHqQ - QVQ 
V(k) 2 



V. 



oo uj — k 
2X1 + W 



-dk 



Hence, 



(E\G(uj)\E) 



uj — k c + in/2 
Auj + in/2 



with 



{Auj - 5 e )(Au + ik j2) - 2\ 2 L ~ A| 

Aw + k/2 
(Auj — uj + )(Auj — ujJ) ' 



Auj — u — k r 



6 e — uj £ k ci 



and 



5 e — ik/2 



± 



\ 



5 e + in/2 



2X1 + XI 



Ft- 



(4.15) 



(4.16) 



(4.17) 

(4.18) 
(4.19) 

(4.20) 
(4.21) 

(4.22) 



In our model, the atoms are in ground state at the beginning and the end. Hence the above 
amplitude is not of direct interest. However, this result is essential for the evaluation of the 
ground-state amplitudes, which is our major objective. The projection of the resolvent by 
Q is 

Q 



QG(uj)Q 



uj - QH Q - QVQ 
9 



u - QH Q - QVQ uj - PH P - PR(uj)P 

V mr, (4.23) 

u - QH Q - QVQ 
yielding the ground-state matrix elements 

<*(*)|<SH|*(V» = ^ + J^^ WE). (4.24, 



These expressions are generally valid 37], and can be applied to any initial photon states. 
In addition, when we are interested in an initial state | \J/) irL satisfying the scattering state 
condition: 

Ve- i6ot \^)- m = for t < 0, (4.25) 

then it can be shown that: 

(Mk)\U{t,0)\Mk')) = 5{k-k')e- ikt [l-2mV{k) 2 (E\G(k)\E)}, (4.26) 

= 5{k-k')e- ikt e i5s{k \ (4.27) 

with 5 s (k) being a fc-dependent real phase shift. Equation (|4.27|) is obtained by noticing that 
the module of the last factor in Eq. 1)4.26)1 is one. While details of the proof will be given in 
the appendix, the physical meaning of Eq. (J4.25)) is obvious. It simply implies that the atoms 
do not experience the field of the injected photon as long as t < 0. Equations 1)4.26)1 and 
()4.27j) become useful when we consider an injected photon prepared in a scattering state. 



V. GENERATION OF ENTANGLEMENT 



After interacting with the atoms, the photon leaks out of the cavity, and its polarization 
can be detected, e.g. by a quarter wave plate and a polarization beam splitter. The atoms 
inside the cavity will then be projected into either the direct product state \LL) or the 
maximally entangled state $ = (\LR) + \RL))/y2, depending, respectively, on whether an 
L or R mode photon is detected. In the following, we consider two different cases where the 
photon in the initial state is (i) prepared in a cavity mode; and (ii) injected from outside 
with a spectral function satisfying Eq. ()4.25j) . 
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A. Cavity photon 

A cavity photon is identified with a quasi-mode photon as defined in Ref. 3^|, which is 
initially created inside the cavity and leaks gradually to the surroundings. In mathematical 
terms, a cavity photon is characterized by the spectral function: 

m = f c (k') = ^ g ;(k% (5.i) 

which corresponds to the cavity line shape. Equations (|5.1j) . (|4.6j) . and (|4.7j) then lead 
directly to 

|*> in = f^dk'^ (V2g L (k')\Mk')} + 9Uk')\Mk')) • (5.2) 

The evolution of the dark states \ip2{k')) is trivial, whereas the resolvent method gives the 
Fourier transform of the time evolution of the states \ipi(k')). As photon detections are 
carried out at the end, the relevant matrix elements of the resolvent are those given by 
Eq. Hence, 

f I (hMfif *\*\ ^ ( 9L{k)f c {k) X L (E\G(uj)\E)V(k) \ 
and from inverse Fourier transform we obtain, as t — ► oo, 

{Mk)m , om » . v-2 (^m + (At _^w_j e -«, ( , 4) 

where all the transients have been neglected, and 

Ak = k - k c . (5.5) 

Notice that although we are looking at the long time behavior, Eqs. ()4.26j) and ()4.27j) do not 
apply to this situation because the initial state given by f c (k') does not satisfy Eq. (|4.25j) . 

The resolvent method only solves the dynamics of the states \ip\{k)). The overall output 
state also includes the contributions of the dark states, in which a trivial phase factor e~ lkt 
is multiplied. With the dark state included again, the output state in the long time limit is 
given by 

ITS lrr\ /*°° J, -iktr A&k - 6 e )(Ak + in/2) - A|,. . 

# out = \LL)® / dke lkt f c {k) y R \k L ) 
J-oo [Ak — uj + )(Ak — u)-) 

n 



where 



a lr = e L -e 



R- 



(5.7) 



If one then detects the photon in the L (R) mode, the atomic state is projected into the 
non-entangled state \LL) (the maximally entangled state with respective probabilities 
Pl and Pr given by: 



fc(k) 



(Ak-S e )(Ak + iK/2) - A| 



(Afc - oo + )(Ak - uj_) 
4A|A|(2A| + XI + k*/2) 



dk 



P, 



R 



2A|A^. 



(2A| + X R ) [(2X1 + A| + k 2 /2) 2 + S^] 

2 

dk 



(5.8) 



fc(k) 



(Ak - uj + )(Ak - lu_) 
_ 4A 2 A|(2A 2 + A| + k 2 /2) 

~ (2A| + A|) [(2A| + A| + ^/2) 2 + « 

If the resonance condition is satisfied, i.e., 5 e = 0, we have 



(5.9) 



R 



4A|A| 



(2A 2 +A|)(2A 2 +A 2 ? + «:V2)- 
In strong coupling regime, where A^, Xr ^> k, 

2AlA/j 



(5.10) 



One can easily prove that 



and equality holds when 



[(2AI + A 2 ,) 



Pr< 



X r = \p2X L . 



(5.11) 



(5.12) 



(5.13) 



Figure El shows the results of Pr versus Xl/Xr when 5 e = 0, and for different cavity 
leakage rates k/Xr = 0, a/2/3, and 7.5. One can see that the probability depends strongly 
on the cavity leakage rate. For finite n, Pr does not attain maximum at Xr = V2Xl exactly 
and the maximum value cannot reach 1/2. In fact, from Eq. ()5.10|) . we have, 



Pr< 



and equality holds when 



X R = V2 



1 + 



2X, 



9 < - 

2 2 



Xr, 



(5.14) 



(5.15) 
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The above results can be better understood by studying the analog of a "three-level 
system" in a perfect cavity. As mentioned previously, our two-atom model reduces to a 
single three-level symmetric system if Eq. (|5.13|) is satisfied. It is well known that such 
a symmetric "three-level system", initially in one of its ground states, oscillates between 
the two ground states. When S e = 0, the oscillations are sinusoidal and complete. When 
there is detuning, the evolution becomes more complicated and perfect population transfer 
from one ground state to the other cannot take place. However, it can be shown that when 
the couplings of the two channels are equal (Eq. (j5.13|) in our model), the time-averaged 
populations of the two ground states are equal, independent of the magnitude of detuning. 
Moreover, as the cavity under consideration is leaky, the probability of detecting an L (or 
R) photon outside the cavity is expected to be proportional to the time the system dwelling 
in each of the ground states. For the case of a good cavity, the upper bound 50% of the 
probability of success is easy to understand, as the system spends equal time on both of the 
ground states in this limit. Hence when the photon leaks outside the cavity, the atom has 



equal probabi 
In Ref. bfi 



ity to decay into an entangled state or a direct product state. 



the hypo-fine levels of a. cesium atom studied in Ret'. [18] were quoted as 



examples for realistic implementation of their scheme. Two different sets of atomic levels are 



studied. In both cases, k/Xr = a/2/3, while Xl/Xr = 1 or y5/6. The numerical results for 
Pr they obtained are 0.43 and 0.45, respectively, which agree with our analytical solution 
shown by the corresponding curve in Fig. El In fact, from Eq. ()5.10j) . our analytic results 
are Pr = 3/7 and 9/20, respectively. 

One of the parameters in Fig. H2 k/Xr = 7.5, is taken from a real CQED setup in 



Ref. |38j. The parameters of the setup, when denoted by the corresponding notations in 
our model, are given by (Xr.tm, ft, 7)/27t = (20, 150, 5)MHz, where 7 is the transverse decay 



rate of the excited state |39j]. (The effect of spontaneous decay has so far been assumed 
negligible in our study. We shall return to discuss this point in Sec. IV Bl ) From Fig. H3 it 
is observed that in this case, the Rabi oscillation frequency is not high enough compared to 
the cavity leakage rate, the bias due to the initial state cannot be erased and this lowers 
the probability of success for poorer cavities. As can be concluded from Eq. (|5.14jl . that in 
the weak coupling limit \l(r) ^ K > Pr — 0, because there is effectively no oscillation at all 
before the photon decays. Hence the scheme using cavity photon is implausible for weak 
coupling cases. However, as will be shown in Sec. IV B[ the entanglement scheme still works 
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nicely in the weak coupling regime with an injected photon. 



B. Injected photon packet 

The motivation to switch from cavity photons to injected photons is two-fold. First, it is 
practically more realistic to consider a photon injected into the cavity than one that already 



exists inside. Second, inspired by the "three-level system" analogy discussed in Sec. IV Al it 
is believed that an injected photon with a duration longer than the cavity leakage time may 
be able to "drive" the system into the entangled state with a higher efficiency We shall 
see that by including general initial photon states other than cavity photon, one may also 
eliminate the requirement for strong coupling. 

Consider an injected photon with a spectral function f(k') satisfying Eq. (|4.25jl . From 
Eqs. (I4.26j) - (j4.27j) . and following similar arguments outlined in Sec. IV Al we have 

/oo />oo 
dke- tkt f{k)C L {k)\k L )-\§)® / dke- lkt f(k)C R (k)\k R } } (5.16) 
-oo J — oo 

where 

(Afc - 5 e )(AP + «V4) - Ak(X R + 2X1) + »«(A£ - 2A|)/2 
L{) (Afc-i«/2)(AJfe-u; + )(AA;-w_) ' { } 

C (k) - V2te^*K\ L \ R 

° Mfcj ~ (Ak - in/2)(Ak - u+)(Ak - u-Y 1 ' 

It is important to note that perfect transfer occurs at the photon frequencies 



Ak = 0, ± V4AI - k 2 /A (5.19) 

because the corresponding amplitudes, 

C L (k) = 0, (5.20) 
C R (k) = e 4A ^, (5.21) 

as long as S e = and Eq. (|5.13J) are satisfied. Therefore, we have an almost unit probability 
of generating the entangled state |$) when f(k) is strongly peaked at one of the frequencies 
defined in Eq. fl5.19|) . The efficiency of entanglement generation depends on the spectral 
width of the injected photon packet, compared with the width of the function C R {k) around 
the peaks. Also notice that detuning only shifts the frequencies in Eq. ()5.19|h and our scheme 
remains plausible by adjusting the peak frequency of the injected photon. 
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We first study an injected photon packet with a complex Lorentzian spectral function 
centered at k c , and with its pole on the upper half plane: 

f(k') = lim - V , ' -e lkT . (5.22) 

Here the parameter Kj n is the spectral width of the input photon packet and r determines 
the initial distance of the photon packet from the cavity. It can be shown that the scattering 
state condition, Eq. ()4.25|) . is satisfied under the limit r — > oo in Eq. (|5.22|) . and hence the 
probability of obtaining the entangled state is given by 

/oo 
\f(k)C R (k)\ 2 dk. (5.23) 
-oo 

Notice that Pr depends only on the module of the spectral function \f(k)\, hence an injected 
photon with spectral function 

f{k') = V ' , (5.24) 

satisfying Eq. (|4.25|) yields the same probability. An injected photon with spectral function 
in the form shown in Eq. 1)5.24)1 can be obtained, for example, by injecting the leaked photon 
emitted by an atom inside another Fabry- Perot cavity. The number of peaks of the spectral 
function depends on the parameters of the cavity and the atom. 

Figure |U shows the dependence of Pr on the spectral width of the injected photon. It is 
shown that in general our method can increase the probability of success to well over 50% 
when the width of the injected photon is sufficiently small compared with the relaxation 
rates of the system. The probability can even approach one in the limiting case of a quasi- 
monochromatic injected photon. 

We point out that in the strong coupling regime, \l(r) ^> the relaxation time of the 
system is Therefore a high probability of success can be achieved when n- m < k. An 
example is shown by the Xl/k = 2.5 curve in Fig. 0J On the other hand, in the bad cavity 
limit, Ai(ij) n, the relaxation rate of the system is modified to (2A| + X r )/k. This 
explains why the Xl/k = 0.13 curve in Fig. 0] requires a smaller Ki n in order to obtain a 
higher probability of success. 

It is useful to note that the curve in Fig. 0] with XjJk = 0.13 corresponds to a realistic 
situation using the parameters in the setup in Ref. j38|. For the cavity photon scheme 
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discussed in the previous section, the probability of success is less than 5%. However, using 
our injected photon scheme, a probability near 50% can still be obtained even if K- m / k ~ 0.3 
(see Fig. HI). 

From Fig. one observes that an intermediate strength of coupling A^/k = 0.25 yields 
best results. This particular value of the parameter is chosen so that the three roots in 
Eq. (j5.19|) coincide. From Eq. (j5.17|) . one sees that if the three roots are far separated, the 
behavior of C*l(A;) around k c is linear: Cx(fc) ~ O(Ak). However, when the three roots 
coincide, the behavior becomes cubic: Ci{k) ~ O(Ak) 3 . As the requirement of energy 
conservation leads directly to |Cx(fc)| 2 + |(7r(A;)| 2 = 1, this means \Cu(k)\ is close to one for 
a wide range of frequencies near Ak = 0. For the specific form of the injected photon we 
assumed, which peaks at the resonance frequency, this means even a larger spectral width 
is allowed to yield a satisfactory efficiency 

The significance of the spectral shape of the injected photon becomes more apparent for 
a Gaussian packet. Consider 



which represents a photon packet with peak frequency k c and spectral width K m . Notice that 



the probability of obtaining the entangled state is again given by Eq. ()5.23|) . 

Figure El shows the dependence of Pr on the width of the injected photon with the same 
set of parameters as in the study of Lorentzian spectral function. One can conclude that 
in general the efficiency is improved using a packet with a Gaussian spectral function. This 
can be understood because for the same spectral width, a Gaussian spectral function is more 
concentrated around the peak frequency than that of a Lorentzian. Remarkably, notice that 
in the case Al/k = 0.25, Pr ~ 1 even for relatively large spectral widths. In fact, we can 
have Pr > 0.999 for /t in ~ 0.3k in a single operation, which can essentially be considered 
deterministic. 

Finally, we would like to address the effects of spontaneous decay. The coupling of the 
atoms to non-cavity modes broadens the excited-levels. This effect can be readily included 
by introducing an imaginary part to the energy of the atomic excited states. In other words, 
we only have to replace u> e byw e — and all the equations remain valid. For cesium 

systems such as that in Ref. |l8j |. we found that the presence of a non-zero 7 only causes a 




(5.25) 



Eq. ()4.25|) is satisfied owing to the phase factor e* fc ' T , and under the limit r — > 00. Hence 
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minute decrease of the probability. For example, by employing the hyperfme levels \L) = 
\QS 1/2 ,F = A,m F = 4), \R) = \6S 1/2 ,F = 4,m F = 2), and |e> = \6P 1/2 ,F = 4,m F = 3), we 
have Xr/Xl = y7/4. With 7 = 0.033/t and X F = 0.25k, one can still obtain a value of Pr 
near 93% using an injected Gaussian photon packet with K m ~ 0.3k. Therefore, the current 
scheme is quite robust against spontaneous atomic decay. 



VI. CONCLUSION 



In conclusion, we study in detail a scheme proposed recently to generate entangled states 
of two identical A-type three-level atoms inside a leaky cavity. The atoms are initially 
prepared in the ground states on the same side of the A systems, with the presence of a single 
photon, which is either in a cavity mode or injected from the exterior of the cavity. First, for 
a cavity photon, we show analytically that an entangled state can be generated in a single 
trial with a certain probability always less than 50%, corroborating the numerical result 
obtained in Ref. [2^. However, their scheme is plausible only in the strong coupling regime. 
By drawing an analogy with a "three-level system" , we provide an intuitive understanding 
of these results. Second, for an injected photon, we show that the probability can be 
increased by injecting a spectrally narrow photon. In particular, we show that an almost unit 
probability can be achieved with a Gaussian packet by exploiting the conditions mentioned 
above. 

Compared with the scheme proposed in Ref. which remains probabilistic for each 
single trial, and requires a feedback mechanism to achieve a unit probability, our proposal 
here can surely increase the success rate in each trial and hence significantly reduces the 
number of trials required to yield an entangled state. We have therefore presented here a 
feasible scheme to generate the maximally entangled state of two A-type three-level atoms 
with a high probability, which can effectively be classified as "deterministic". This is also 
in contrast to recent probabilistic schemes, in which intrinsic uncertainties are inherent 



in quantum measurement processes 



24L |26| . However, the main challenge of realizing our 



scheme is the requirement of a single photon source, which is now under active investigations 
and may become feasible in the near future j^j. 
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APPENDIX A 



In our model, we are interested in initial photon states that may or may not be obtained 
by using Eqs. ()4.26|) . For example, it can be shown that with a cavity photon defined by 
Eq. (|5.1|) . or a Gaussian packet defined by 

2(k' - k c ) 



f(k') = i?r ^ — exp 



(Al) 



Eqs. (I4.24|) and ()4.26|) lead to different results. However, by introducing a phase factor in 
the spectral functions, Eq. (J4.26|) becomes valid for a Gaussian packet defined by Eq. ()5.25|) 
or a "displaced" cavity photon in Eq. (|5.22jl . In addition, for some packets without a phase 
factor like e lk ' T , such as a complex Lorentzian defined by Eq. (J5.24J) . Eq. ()4.26j) still yields 
the correct results. Hence, it is desirable to obtain a necessary and sufficient condition for 
the applicability of Eq. ()4.26|) . Here we will prove that Eq. ()4.25|) is the condition we are 
looking for. 

For the initial state 

/oo 
dk'f(k')\Mk% (A2) 
-oo 

Eq. (l4~2"oTl leads to 

IZo f{k')V{k')e~ ik,t dk' = for t < 0. (A3) 

Equation (|A3J) implies 

/OO roo roo 

dk'f{k')V{k')e iiz - k)t = dt dk'f(k')V{k')e i{ - z ~ k)t . (A4) 
-oo JO J — oo 

The proper meaning of the above equation must be defined by certain limiting process. For 
our problem, we choose to introduce to the real number z a small imaginary part, which has 
the same sign as t. Hence 

roo roo roo roo 

lim / dt dk'f(k')V(k')e {iz - ik ~ e)t = lim / dt dk'f{k')V{k')e {iz - ik - sgIimt . 

e^0+ Jo J-oo e^0+ J -oo J -oo 

(A5) 
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Therefore 

/oo j yOG 

f(k f )V(k')—^— — -dk' = 2nj f(k')V(k')S(z-k')dk', (A6) 

or 

lim f° f{k ' )V{k l dk' = -2nif(z)V(z). (A7) 

e^0+ J-oo z + le — k' 

From Eq. ([4.24)1 . we can conclude that 

/CO . r /CM*' 



t+ M WCM | {| rwzi it - (A9) 

lu — k to — k J-oc to — k' 

For t > 0, the evolution is governed by the retarded propagator. We can thus let uj approach 
the real axis from above, and from Eq. (|A7[) . we have 



W1)fe |,),4- 2 „«»W. (AID) 

Because f(k) represents a wave-packet, V has finite support in space, and (E\G(u)\E) 
governs the decaying process of the excited state, k is the only real pole that one needs to 
keep in the asymptotic long time limit. Hence by performing the inverse Fourier transform, 
we have 

(Mk)\U(t,0)\V) in = f(k)e- m (l - 2mV(k) 2 (E\G(k)\E)) , (All) 
which is equivalent to defining the S'-matrix as shown in Eq. (|4.2fij) . 
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FIG. 1: A sketch of the system in our model. Two A- type three- level atoms are trapped near 
the center of a one-sided leaky cavity. The continuous frequency modes shown in the figure are 
adopted to study the generation of atomic entangled state by a cavity mode photon or an injected 
photon packet. 
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FIG. 2: A schematic plot of |t/fc(xo)| 2 versus the frequency A; at a point xq inside the cavity. Each 
peak corresponds to a quasi-mode of the cavity. The separations and widths of the quasi-modes 
are given by it /I and k, respectively. In our model, we assume that the atomic frequency is close 
to one of the quasi-mode frequencies, and adopt the single mode approximation. 




FIG. 3: A plot of probability of successful generation of the entangled state, Pr, versus Xl/^r as 
obtained from Eq. (|5.9j) . with 5 e = and k/\r = 0, y/2/3, and 7.5. 
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FIG. 4: A plot of probability of successful generation of the entangled state, Pr, versus K- m / 'n for an 
injected photon with a Lorentzian spectral function given by Eq. (|5.24|) . with 5 e = 0, = 2.5, 

0.25, and 0.13. 
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FIG. 5: A plot of probability of successful generation of the entangled state, Pr, versus k^/k for 
an injected photon with a Gaussian spectral function given by Eq. (|5.25|) . with 5 e = 0, Xl/k = 2.5, 
0.25 and 0.13. 
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